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Abstract. The symmetric Anderson impurity model with a hybridization vanishing 
at the Fermi level, Aj oc \oj\ t , is studied via the numerical renormalization group 
(NRG) at T — 0; and detailed comparison made with predictions arising from the 
local moment approach (LMA), a recently developed many-body theory which is 
found to provide a remarkably successful description of the problem. Results for 
the 'normal' [r = 0) impurity model are obtained as a specific case, and likewise 
compared. Particular emphasis is given both to single-particle excitation dynamics, 
and to the transition between the strong coupling (SC) and local moment (LM) phases 
of the model. Scaling characteristics and asymptotic behaviour of the SC/LM phase 
boundaries are considered. Single-particle spectra D(lu) are investigated in some detail, 
for the SC phase in particular. Here, in accordance with a recently established result, 
the modified spectral functions T(uj) oc |a;| r Z?(w) are found to contain a generalized 
Kondo resonance that is ubiquitously pinned at the Fermi level; and which exhibits 
a characteristic low-energy Kondo scale, wkW, that narrows progressively upon 
approach to the SC— >LM transition, where it vanishes. Universal scaling of the spectra 
as the transition is approached thus results. The scaling spectrum characteristic of the 
normal Anderson model is recovered as a particular case, that exemplifies behaviour 
characteristic of the SC phase generally, and which is captured quantitatively by the 
LMA. In all cases the r-dependent scaling spectra are found to possess characteristic 
low-energy asymptotics, but to be dominated by generalized Doniach-Sunjic tails, in 
agreement with LMA predictions. 
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1. Introduction 

The Anderson impurity model [1], in which a correlated spin-| impurity is coupled to, 
and quenched by, the low-energy excitations of a non-interacting metallic host, has long 
occupied a central role in condensed matter theory (for a comprehensive review see e.g. 
[2]). Its low-energy properties — those of a conventional Fermi liquid — are of course 
contingent upon coupling to a metallic host, whose density of states phost(^) is non- 
vanishing at the Fermi level, uj = 0. However the 'normal' impurity model is itself a 
limiting case (r = 0) of a more general class, in which the host spectrum contains a 
soft-gap at the Fermi level, phost(^) oc |c<j| r with r > 0. Such models, exhibiting a much 
wider spectrum of physical behaviour than the r = limit (as discussed below), are 
naturally of considerable theoretical interest; but not solely so, there being a number 
of potential candidates for soft-gap behaviour, including various semiconductors [3] and 
heavy Fermion superconductors [4], certain two-dimensional systems [5,6] and quasi- 
one-dimensional metals described by a Luttinger model [7]. 

The soft-gap problem was first considered by Withoff and Fradkin [8] in relation 
to the corresponding Kondo model. Much study of both soft-gap Kondo and Anderson 
models has since ensued, in particular via poor man's scaling [8,9], large- N expansions 
(with N the impurity degeneracy) [8,10,11], the numerical renormalization group (NRG) 
[9,12-15] and perturbation theory in the impurity interaction strength, U [16]. From 
such, it is now well established that two distinct types of ground state arise, between 
which in general a quantum phase transition occurs at a critical, finite host-impurity 
coupling: (i) a doubly degenerate local moment (LM) state in which the impurity 
spin remains unquenched; and (ii) a strong coupling (SC) state in which the impurity 
spin is locally quenched and a Kondo effect manifest, and whose properties have been 
argued to represent a nontrivial but natural generalization of the Fermi liquid behaviour 
ubiquitous to the normal Anderson model [15,16]. The underlying physics is known to 
be particularly rich for the particle-hole symmetric problem, to which NRG studies 
in particular have devoted considerable attention, encompassing both thermodynamic 
properties [12-15] and impurity single-particle spectra [14]. 

For impurity models generally, a theoretical description of dynamics, and in 
particular single-particle dynamics, poses a significant challenge. Even for the normal 
spin-| Anderson model, current theories such as the non-crossing approximation [17- 
20], 1/N expansions [21-23] or slave boson approaches [24-26] are well known to possess 
significant limitations; whether it be an inability in practice to handle finite interaction 
strengths, to describe excitation dynamics on all energy scales, or to recover low-energy 
spectral characteristics (Fermi liquid behaviour, for instance, or Doniach-Sunjic tails [27] 
in the 'wings' of the Kondo resonance). Benchmark results are of course now accessible 
via NRG calculations (see e.g. [14,28,29]), or by a combination of Quantum Monte Carlo 
calculations and the maximum-entropy method (at finite temperature) [30] . It is against 
such that the accuracy, qualitative or otherwise, of extant theories — or new ones - 
must be assessed. 
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In this paper we consider the symmetric soft-gap Anderson model, including the 
'normal' limit of r = 0. Our primary aim is to make detailed comparison between NRG 
calculations and a new theoretical approach to the problem, known as the local moment 
approach (LMA), that has recently been developed to treat the soft-gap Anderson model 
in [31] (denoted hereafter as I); and from which results for the normal model, to which 
the approach was originally applied [32], are recovered smoothly in the limit r — > 0. 

The dominant focus of the LMA is a treatment of single-particle dynamics, 
embodied in the impurity single-particle spectrum D{uj) — on all energy scales and 
for any interaction strength — although an essential element of it also permits the 
SC/LM transition, and associated phase boundaries, to be addressed directly. The 
LMA is naturally non-perturbative, with both the notion of local moments and the 
possibility of either a SC or LM state introduced explicitly and self- consistently from 
the outset, as reflected in the adoption of an underlying two-self-energy description for 
the impurity Green function. The approach itself is introduced briefly in §2 (full details 
being given in I); where in addition we emphasize the importance of a [/-independent 
pinning theorem for the single-particle spectrum of the SC phase at the Fermi level, 
recently established on general grounds in [16]. We then outline a number of predictions 
arising from the LMA that, in addition to single-particle spectra, can be tested directly 
against NRG calculations. These include two regimes of one-parameter scaling for the 
SC/LM phase boundaries; asymptotically exact results for the small-r behaviour of both 
the phase boundaries and Kondo scale characteristic of the SC phase; and the prediction 
that, as for the normal r = Anderson model, the so-called modified spectral functions 
J-{oj) oc \u\ r D{uj) should exhibit scaling behaviour upon approach to the SC— ^LM phase 
boundary — where the Kondo scale vanishes — resulting thereby in an r-dependent 
family of universal scaling spectra. 

The NRG procedure employed is outlined in §3, with particular emphasis on how 
information on ground state properties can be extracted, and how they depend upon the 
parameters of the NRG calculations. In §'s 4-6 we turn to detailed comparison between 
NRG and LMA results. Phase boundaries, including their scaling characteristics and 
asymptotic behaviour, are considered in §4; comparison with previous NRG calculations 
[14,15] is also made. Single-particle spectra are considered in §5, beginning with direct 
comparison between NRG and LMA results for the 'bare' D(u), for both SC and LM 
phases. This is followed by consideration of the modified spectra ^{to) , which have been 
argued (see I and [16]) to provide a much more revealing description of single-particle 
dynamics characteristic of the SC phase: in direct parallel to the normal Anderson 
model, the J-'iw) — for which the pinning theorem [16] is confirmed numerically — are 
seen to contain a generalized Kondo resonance whose width is proportional to the Kondo 
scale wk = wkWj progressively narrows as the SC^LM transition is approached, and 
vanishes at the transition itself. 

The latter leads naturally to the issue of spectral scaling, to which we turn in §6. 
We begin by considering the normal Anderson model, for which the strong coupling 
(large-f/) NRG scaling spectrum D = D(lu/lu k ) is first obtained, and compared to the 
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LMA result [32]. The quantitative agreement between the two is quite remarkable; and 
in qualitative terms, while characteristic Fermi liquid behaviour is naturally recovered 
for |c<j|/c(Jk << 1, the scaling spectrum is clearly dominated by the Doniach-Sunjic 
(DS) power-law tails [27,28,30] indicative of an incipient orthogonality catastrophe. 
For r > 0, universal scaling of the SC J-{uS) is likewise found from NRG calculations. 
The r-dependent scaling spectra are extracted, and indeed found to possess both the 
characteristic low- frequency (|cj|/c<jk << 1) asymptotics and generalized DS tails that 
the LMA predicts. 

Our overall conclusion from this study is that the LMA provides a rather successful 
description of the soft-gap and normal Anderson models, and appears to transcend many 
limitations of extant theoretical approaches. Its specific predictions are well borne out 
by benchmark NRG calculations; the level of agreement with which, both qualitative 
and quantitative, provides encouraging impetus to further study and development of 
the basic approach. 



2. Background 

In standard notation, and with the Fermi level taken as the origin of energy, the 
Hamiltonian for the spin-| Anderson model is given by 

H -^host -^impurity -^hybridization 

= ^2 tkflka + + \ U hi-a)n™ + ^2 V ik(4a C ka + c{ a C ia ) (2.1) 

k,a a k,cr 

where e& is the host dispersion (with corresponding spectrum Phost(^)), Mfc is the 
hybridization matrix element and e x the impurity level; for the symmetric model 
considered here, q = — £7/2 with U the on-site Coulomb interaction. The host-impurity 
coupling is embodied succinctly in the hybridization function A (a;) = J^ fc |V;fc| 2 /(cj + 
i?7Sgn(u;) — e*.) (where r\ = 0+), such that A (a;) = Ar(cj) — isgn(a;)Ai(a;) with 

A I (uj) = nJ2\y^\ 2 S(io-e k ). (2.2) 
fc 

For the soft-gap model, Ai(u) takes the power-law form 

Ai(u) = A \uj\ r 6(D - M) (2.3) 

with r > and bandwidth D (6(x) being the unit step function). Note from Eq. 

(2.2) that the form Eq. (2.3) arises from a constant V\k and a soft-gap host spectrum 

Phost(^) — Ylk^( u ~ e fc)' although a separate specification of the {Vik} and host 

eigenvalues {e^} is not in fact required to specify Ai(u>). The corresponding real 

part A r (cj) follows by Hilbert transformation, its low-u; behaviour being given by 

Ar(w) = — sgn(c;)tan( : |r)Ao|a;| r + 0{\u\/D) for r < 1 (as here considered in practice). 

The problem is thus characterized by three bare energy scales, viz U, D and (as 

i 

convenient) either AQ~ r or A D r ; one of which may be used as the energy unit, resulting 
in two independent dimensionless parameters. 
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The normal Anderson model, r = is recovered as a particular limit of Eq. (2.3). 
Here, quenching of the impurity spin is ubiquitous and the system well known (see 
e.g. [2]) to be a conventional Fermi liquid for all U > and any Ao > (Ao = 
corresponding to the atomic limit where the impurity/host trivially decouple). For 
r > by contrast, two distinct possible ground states are known to exist [9,14,15] — a 
doubly degenerate local moment (LM) state, and a strong coupling (SC) state in which 
the impurity spin is locally quenched and a Kondo effect manifest — between which 
in general a nontrivial quantum phase transition occurs at a critical coupling strength. 
This will be investigated in the following sections, employing both the NRG and LMA. 
Details of the former will be given in §3; first, we outline the local moment approach 
and its predictions [31,32] relevant to the present work. 

2.1. Dynamics and the LMA 

The LMA focuses explicitly on single-particle dynamics, embodied in the zero 
temperature impurity Green function G(uj) («-> G(t) = — i(T{ci a (t)c ia })) and hence 
single-particle spectrum D{uj) = — 7r -1 sgn(u;)ImGr(ti;); and G{uS) is usually expressed as 

G(u) = [lu + ir?sgn(ccO - A(u) - S(cj)] -1 (2.4) 

which simply defines the conventional single self-energy £(c<->) [33]. As far as dynamics are 
concerned, our primary interest in this paper resides in the SC phase with its attendant 
low-energy Kondo scale. In this regard we first note an important conservation upon 
the single-particle spectrum at the Fermi level, ou — 0, hitherto established in [16] on 
general grounds: namely, defining a modified spectral function J-{u) by 

F(u) = ttA [1 + tan 2 (f r)] \u\ r D(u), (2.5) 

that 

F( u = 0) = 1 (2.6) 

throughout the SC phase. For r = 0, this recovers as a particular case a well known 
result for the normal impurity model, usually viewed as a consequence of the Friedel 
sum rule (see e.g. [2]): that ttAqD(uj = 0) = 1 — the spectrum is always pinned at 
the Fermi level uo = 0. Equation (2.6) generalizes this pinning condition to arbitrary r 
for a SC state and, as discussed in [16], embodies the fact that many-body interactions 
have no influence in renormalizing the low-u; behaviour of D{uj) in the SC phase: its 
uj — > asymptotic form D(u>) ~ M~ r is precisely that of the non- interacting limit, 
consistent with the view [15,16,31] that the r > SC state constitutes a non-trivial 
but natural generalization of Fermi liquid physics. The extent to which the pinning Eq. 
(2.6) is captured in practice should also provide a good test of the accuracy of NRG 
calculations, as considered in §5. 

In seeking to describe single-particle dynamics for both SC and LM phases, the 
LMA [31,32] eschews direct calculation of the single self-energy T,(u), and instead 
employs a two- self-energy description with G(u) expressed formally as 

G(u) = ±[Gi(u) + G i (u)] (2.7a) 
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where 

G a {uj) =L + i V sgn{uj) - A(u) - * (2.7b) 

(and a =| / | or +/—). The interaction self-energies So-(c<j) (= — — uj) by particle- 
hole symmetry) are separated as [33] 

E <T H = -f?7H + E <T H (2.8) 

into a purely static Fock contribution (with local moment |/i|) that alone would be 
retained at the simple mean-field level of unrestricted Hartree Fock (UHF), together with 
an u;-dependent contribution T, a (ou) containing the dynamics that, at low frequencies 
in particular, are naturally central to the problem. There are several reasons for 
adopting a two- self-energy description. First, it is a physically natural description of the 
doubly degenerate LM state, the self-consistent possibility of which must be introduced 
explicitly from the outset if one seeks to devise a non-perturbative approach that can 
simultaneously handle the possibility of both LM and SC states. Second, it provides a 
tangible means of developing a many-body approach to the problem that starts from, 
but successfully transcends the limitations of, the static mean-field description (UHF). 
Finally, conventional perturbation theory (PT) in U about the non- interacting limit, 
which underpins traditional diagrammatic approaches to the single self-energy £(u;), is 
at best limited: for example the existence of a SC/LM transition at a finite interaction 
strength for < r < \ [9,14,15] attests to a finite radius of convergence for conventional 
PT, while for | < r < 1 there is evidence to suggest a vanishing radius of convergence 
for PT in U [16]. 

The LMA, described in detail in I and [32], has in practice two essential elements. 

(i) It includes in the self-energies S ct (cj) a non-perturbative class of diagrams (Fig. 3 of I) 
that embody dynamical coupling of single-particle excitations to low-energy transverse 
spin fluctuations, and hence capture the spin-flip scattering essential to describe the 
Kondo, or spin-fluctuation, regime. Other classes of diagrams, involving primarily 
charge and longitudinal spin fluctuations, may also be included in the S ct (cj) (see e.g. 
Fig. 9 of [32]); but retention of the dynamical spin- flip scattering processes is essential. 

(ii) In describing the SC phase (and hence delineating the boundaries thereof), the 
spectral pinning Eq. (2.6) at the Fermi level u = precisely, is enforced as a self- 
consistency condition; a necessary/sufficient condition for which is readily shown (§5.1 
of I) to be that (Ret a (co = 0) =) = 0) = 0, i.e. £*(cj = 0) = § U\/i\ from Eq. 
(2.8). But S ct (cj) depends both explicitly on U and parametrically on (since the 
bare propagators that enter Yi a (u)) are of mean-field form and thus depend upon the 
static mean-field self-energy — Hence, enforcement of the u = spectral pinning 
characteristic of the SC state is guaranteed by 

£> = 0;C/;f[/H) = f[/|/x|. (2.9) 

For a chosen r and U, Eq. (2.9) is a self-consistency equation for the local moment 
the limits of stability of solutions to which, upon increasing U, thus yield the 
SC/LM phase boundary (which is thereby found, correctly, to be coterminus with that 
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obtained upon approach from the LM phase, see §6 of I). Most importantly, as detailed 
in I, self-consistent solution of Eq. (2.9) introduces naturally into the problem a low 
energy scale that (obviously) has no counterpart at crude mean-field level: the Kondo, 
or spin-flip scale uj m = u; m (r), whose physical significance within the LMA is twofold. 
First, it corresponds to the energy cost to flip the impurity spin (as manifest in a strong 
resonance centred upon u = u m in the transverse spin polarization propagator [31,32]). 
With increasing interaction strength in the SC phase, u m (r) diminishes progressively 
and vanishes as the SC/LM phase boundary is approached; while w m (r) = throughout 
the LM phase — as expected physically, since for a doubly degenerate LM state with 
finite weight on the impurity there is no energy cost to flip a spin. Second, in the SC 
phase, cu m (r) > sets the finite timescale (~ l/a> m ) for restoration of the locally broken 
symmetry inherent to the zeroth-order mean-field level of description, as manifest for 
example in the fact that for u/u m <C 1, Sf(o;) and Xj,(c<-0 coincide with each other (and 
hence with the conventional single self-energy E(u;)); while for the LM phase by contrast, 
where u m = 0, there is naturally no such symmetry restoration and S|(cj) ^ Ej(^) even 
as lu — > 0. 

The preceeding comments are merely intended to provide a brief overview of the 
strategy behind the LMA, and full details are given in I. We now highlight some 
predictions of the approach that will be tested against NRG calculations in §'s 4-6, 
beginning with the SC/LM phase boundaries. 

In qualitative terms, and for any U > 0, the LMA yields the existence solely of 

LM states for all r > \. This agrees with extant NRG calculations [14,15] (although 

we note that for U = precisely the ground state is known [16] to be SC for all 

< r < 1; so for \ < r < 1 there is a SC/LM transition 'at' U — itself, a 

fact that underlies the breakdown of conventional PT in U in this r-regime, see [16] 

and I). Moreover, for the normal Anderson model r = 0, the LMA correctly yields 

a normal Fermi liquid for all U > and A > [32] - a well known fact, but 

one by no means guaranteed a priori in an approximate theory. The SC/LM phase 

boundary is thus confined in effect to < r < |, and its r-dependence will in general 

depend upon the two independent dimensionless parameters that, as mentioned above, 

i 

characterize the model: one of A D r /U or A^ r /U, together with U/D. However the 
LMA predicts two distinct regimes of one-parameter scaling of the phase boundaries, 
according to whether [//2 > D or < D — corresponding physically to an impurity 
level \ei\ = |C/ that lies respectively well outside or within the host conduction band 
of width D: (a) For [//2 > D, the phase boundary depends universally on A D r /U 

alone, and not upon U/D; while (b) for (7/2 « D the transition line depends solely 

i 

on A U r /U = (A 1 ~ r /£/) 1 ~ r . The latter regime is known to exist from previous NRG 
studies [15], the former will be investigated via the NRG in §4. 

The small r behaviour of both the phase boundary, and the vanishing of the Kondo 
scale u> m (r) as the transition is approached from the SC phase, are naturally important. 
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These are given respectively within the LMA by (§6.1 of I) 

A X r \ r ^0 7T 



U l c ~ 8 r {2M) 



and 



i- 17 



Ur)\ (2 ' U) 
where A = min[D, C//2]. We believe the exponent of 1/r, and Eq. (2.10) for the phase 
boundary, to be asymptotically exact as r — > 0. The latter is seen by taking the 
limit r -> of Eq. (2.11), using Eq. (2.10), to yield u m {r = 0) = Aexp(-7rC//8A ). 
This is the strong coupling (large-C/) Kondo scale for the normal Anderson model. 
The prefactor of A is merely an approximate high-energy cut-off, but the exponent 
of — 7r£//8A is exact, as known from the Bethe ansatz solution [34] and poor man's 
scaling (see e.g. [2]); and recovery of it hinges on the asymptotic validity of Eq. (2.10) as 
r — > 0. Equations (2.10,11) will likewise be compared to NRG calculations in §4. Note 
further from Eq. (2.10) that the condition for mapping the symmetric Anderson model 
onto the corresponding Kondo model under a Schrieffer- Wolff transformation [35], viz 
A A r /U <C 1, is evidently satisfied along the phase boundary line as r — > 0; using Eq. 
(2.10) the critical J c for the Kondo model as r — > can thus be obtained and, as detailed 
in I, recovers precisely the result obtained by Withoff and Fradkin [8] from poor man's 
scaling. 

One strength of the LMA is its ability to describe single-particle dynamics, which 
are compared to NRG calculations in §'s 5 and 6. The known low-cu spectral signatures 
of the LM and SC phases, viz D(oo) ~ \oo\ r and \oo\~ r respectively as uo — > [14], are 
correctly recovered by the approach; and it gives excellent agreement with extant NRG 
results for D{uS) [14] (see Fig. 15 of I). However, as shown in [16] and I, the most revealing 
expose of single-particle dynamics in the SC phase is evident not in D(u) itself, but in 
the modified spectral function J-{u) oc \uo\ r D{uo) (Eq. (2.5)) wherein the unrenormalized 
\uo\~ r divergence in D(ui) has been removed. In thus 'exposing' the low-frequency many- 
body renormalization characteristic of the Kondo effect, the Kondo scale is now directly 
manifest in J-(u): in addition to spectral pinning at the Fermi level u = (Eq. (2.6)), 
J-{oj) contains a generalized Kondo resonance with a width proportional to cu m (r), that 
narrows upon progressive increase of U in the SC phase and vanishes at the SC/LM 
phase boundary; see §8 of I. This behaviour, characteristic of the SC phase for all r < |, 
is in obvious parallel to that for the normal Anderson model, and is further support for 
the notion of the SC phase as a 'generalized Fermi liquid'. NRG results for J-{uj) will 
be given in §5. 

Most significantly, since the width of the Kondo resonance in ^{oo) vanishes as the 
SC— >LM transition is approached (U — > U c (r)—), the above comments suggest — and 
the LMA indeed predicts (see §8 of I) — that J-{uS) should exhibit scaling behaviour 
as the phase boundary is approached from the SC side: i.e. T(uj) should be a universal 
function of uj/u) m (j-) for any r where a SC state exists. An r-dependent family of 
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universal scaling spectra is thus expected to arise, encompassing as a specific case the 
well known scaling for the normal Anderson model r = as U — > oo (= U c (r — 0), see 
Eq. (2.10)). This is considered in §6 where we first compare the NRG and LMA scaling 
spectra for the normal Anderson model. Further, the LMA predicts the universal .F(u;)'s 
to exhibit characteristic r-dependent low-frequency behaviour (for u/uj ni < 1), as well 
as Doniach-Sunjic tails [27] that are a reflection of the orthogonality catastrophe and 
represent a generalization of those known to arise in the r = Anderson model [28,30]; 
this too is compared to NRG calculations in §6. 

3. NRG approach 

Detailed discussion of how the NRG can be applied to the soft-gap Anderson model can 
be found in [9,12-15]. Here we mention briefly only those aspects of the approach 
necessary to understand how e.g. information on ground state properties can be 
extracted, and how they depend upon the parameters of the NRG calculations. 

The NRG is based quite generally on a logarithmic discretization of the energy axis, 
i.e. one introduces a parameter A and divides the energy axis into intervals [A~( n+1 ), A _n ] 
for n — 0, 1, oo [36,37]. With some further manipulations [36,37] the original model 
may be mapped onto a semi-infinite chain, which can then be solved iteratively by 
starting from the impurity and successively adding chain sites. The coupling between 
two adjacent such sites n and n + 1 vanishes as A _ri//2 for large n, whence the low-energy 
states of the chain with n + 1 sites are generally determined by a comparatively small 
number N s of states close to the ground state of the n-site system. In practice one 
retains only these N s states from the n-site chain to set up the Hilbert space for the 
n + 1 site chain, thus preventing the usual exponential growth of the Hilbert space as n 
increases. Eventually, after unrg sites have been included in the calculation, addition of 
another site will not change significantly the spectrum of many-particle excitations; the 
spectrum is very close to that of the fixed point, and the calculation may be terminated. 

An example of such a flow diagram for some of the lowest lying energy levels as 
a function of the chain length, n, is given in Fig. 1 for the soft-gap Anderson model 
with r = 0.4, U/D = 10~ 3 and A D r /D = 0.0075; the parameters used for the NRG 
calculations are A = 2 and N s = 300. The states are labelled by the quantum numbers 
Q (which characterizes the number of particles measured relative to one particle per 
site), and the total spin, S. As mentioned above, the energy scale is reduced in each 
step by a factor A 1 / 2 . To allow for a direct comparison of the energies for different chain 
lengths, it is thus convenient to plot A n ^ 2 E n instead of the eigenvalues E n of the n-site 
chain directly. As is apparent from Fig. 1, the properties of the system in this case do 
not change further for chain lengths n^Rc > 120. Without going into details here, one 
can show that the distribution of energy levels for n > 120 in Fig. 1 is characteristic of 
the SC phase of the model. 

If by contrast we choose instead a value of A D r / D = 0.006, we obtain the flow 
diagram shown in Fig. 2. Here it is evident that the fixed point level structure is entirely 
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Figure 1. NRG flow diagram for the lowest lying energy levels for AqD t / D = 0.0075. 
The system flows to the SC fixed point. 
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Figure 2. Flow diagram for the lowest lying energy levels for AqD t /D = 0.006. 
The system flows to the LM fixed point. 

different from the SC solution, and indeed this particular pattern is now characteristic 
of the LM phase of the model. We can thus conclude, simply from inspection of the 
two flow diagrams, that the critical (AqD t / D) c separating the SC and LM phases of 
the soft-gap Anderson model for the model parameters specified, lies in the interval 
[0.006,0.0075]. Performing a whole series of such NRG calculations leads to Fig. 3 
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(illustrated for the first excited state with Q=l and S=0), and the conclusion that 
the critical A D r /D lies between 0.00670 and 0.00675. This method thus enables in 
principle a determination of the critical coupling to arbitrary precision. Note further 
from Fig. 3 that, in the vicinity of the critical coupling, the approach to the SC or 
LM value of A n / 2 E n occurs at progressively higher values of n. This corresponds to 
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Figure 3. Flow diagram for the first excited state with quantum numbers Q = 1 
and 5 = 0, for various values of A = AqD t /D. 
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Figure 4. A-dependence of the critical A' c = (AqD v / D) c for U/D =0.1 and 
r = 0.1 (circles) and r = 0.4 (squares). The A' c for r = 0.1 is multiplied by a 
factor of 10. 



the vanishing of a low-energy scale as the transition is approached from either side, the 
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low-energy scale being proportional to the Kondo scale when coming from the SC phase. 
A low-energy scale may be defined as cu* = Ar n *l 2 with n* given by \x n<t — x^] = 0.01, 
with x n = K a l 2 E n . In the SC phase the low-energy scale can also be determined via e.g. 
the width o;k of the Kondo resonance (see §5), and the NRG results show that uu* oc c^k 
(the two are not of course strictly equal, since a low-energy scale is determined only up 
to a pref actor). 

It is of course obvious that for any A > 1 the NRG constitutes an approximation 
to the continuum system, but becomes exact in the limit A — > 1. Performing this limit 
directly is naturally not possible since one must simultaneously increase the number of 
retained states, N s , to infinity; but one can study the A and iV s -dependences of the NRG 
results and perform the limits A — > 1, N s — > oo by extrapolation. As an example, Fig. 
4 shows results for the A-dependence of the critical (AqD t / D) c for a fixed N s = 300; 
for both r = 0.1 and r = 0.4, and with U/D = 0.1. It is evident from Fig. 4 that 
the A-dependence is rather mild, although increasing slightly as r — > |. Similar study 
of the A^-dependence also reveals a rather weak variation that again does not depend 
significantly on r. From such considerations we conclude that for practical purposes, a 
choice of A = 2 and N s = 300 will in general be sufficient to determine e.g. the SC/LM 
phase boundaries to good accuracy and with minimal numerical effort. One should 
nonetheless keep in mind that due to the A-dependence the actual critical coupling 
(A D r I D) c will always be a few percent smaller, the deviation increasing mildly as 
r -> \. 

Finally, we remind the reader that the NRG also enables calculation of dynamical 
properties such as the single-particle spectrum D(u), as detailed in [14] and [38]; and 
results for which will be given in §'s 5 and 6. 

4. Phase boundaries 

In discussing phase boundaries between SC and LM states, we consider first the case 
U/D <C 1 that has been the focus of previous NRG studies [14,15]. Here the phase 
boundary for given r depends solely on the ratio A U r /U, a result known to arise from 
application of poor man's scaling to the soft-gap Anderson model [9]. It may however 
be understood quite simply by noting that U/D — > corresponds to the conventional 
limit of an infinite host bandwidth, D — > oo (which may be taken with impunity 

for r < 1, see e.g. I); in which limit the D-scale drops out of the problem, and the 

i 

model thus depends solely on the dimensionless ratio U = U / A^ r . Phase boundaries 
correspond to a critical U c (r), an alternative but equivalent way of displaying which is 
via (A U r /U) c = Ur\r). 

The above behaviour for U/D <C 1 has been observed in NRG calculations by 
Gonzalez-Buxton and Ingersent (G-BI) [15], who find in practice excellent scaling of the 
phase boundary for U/D = 0.2 and 0.02 (see Fig. 20 of [15], where the quantity p$J c 
there plotted for the symmetric model is (8/7r)(AoA r /£/) c with A = U/2). As noted in 
I however, the phase boundary (A U r /U) c obtained originally by Bulla, Pruschke and 
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Figure 5. For U/D <C 1, comparison of NRG phase boundary (AoU r /U) c versus r 
(solid diamonds) with that obtained via the LMA (solid line). Previous NRG results 
[15] are also shown (open squares), together with the r — > Kondo asymptote of nr/8 
(dashed line). 



Hewson [14] for U/D = 10~ 3 differs significantly from that found by G-BI (see Fig. 9 

of I for an explicit comparison). We have thus repeated a determination of the phase 

boundary for U/D = 10~ 3 using the more refined procedure outlined in §3. As shown in 

Fig. 5, good agreement with the G-BI results is now found: the difference between the 

two sets of data is small, and readily attributable to the choice of A = 2 and N s = 300 

in the present calculations, which (see §3) will yield a slight overestimate of (A U r /U) c 

(the G-BI results [15] having already been extrapolated to the continuum limit). 

Fig. 5 will be discussed further below, before which we consider the opposite 

limit of U/D ^> 1. Here, in contrast to U/D <C 1, the LMA predicts (see §6.3 of 

i 

I) that none of the bare scales U, D or A^ r drop out of the problem, but that the 

phase boundary exhibits one-parameter scaling when expressed in terms of AoD r /U(= 
i 

(U/A^ r ) r ~ 1 (D/U) r ). This too is physically understandable, since for \e{\ — U/2 D 
where the impurity level lies well outside the host band, the impurity-host coupling is 
controlled by the hybridization A\(D) = AqD t which, together with U, sets the natural 
energy scales upon which ratio the phase boundary thus depends. This is confirmed 
by NRG calculations, as shown in Fig. 6. In Fig. 6a we show the critical (A U r /U) c 
versus r for U/D = 10 3 , 10 2 , 10 and (as in Fig. 5) 10~ 3 ; expressed in this form, the phase 
boundary curves depend strongly upon the chosen U/D. For the same U/D ratios, Fig. 
6b by contrast shows the critical (AqD t /U) c versus r, from which the predicted one- 
parameter scaling for U/D >> 1 is seen clearly; in particular, the results for U/D = 10 2 
and 10 3 are essentially indistinguishable. 

In Fig. 7, NRG results for (A D r /U) c versus r (with U/D = 10 3 ) are compared to 




Figure 6. NRG phase boundaries for U/D = 10 3 (squares), 10 2 (circles), 10 (triangles) 
and 10 -3 (diamonds). Expressed as (a) (A U r /U) c the curves are strongly [//in- 
dependent; but as (b) (Ao-D r /[/)c, the predicted one-parameter scaling is evident for 
U/D 3> 1. Lines connecting points are a guide to the eye only. 



those arising from the LMA for U/D ^> 1; while for U/D <C 1, Fig. 5 likewise shows 
the corresponding LMA result for (AoD r /U) c . In both Figs. 5 and 7 we also indicate 
the predicted r — > asymptotic behaviour of the phase boundary that is symptomatic 
of the Kondo limit, viz nr/8 (Eq. (2.10); behaviour that is indeed recovered from the 
NRG calculations, regardless of whether [7/D > 1 or < 1. In the latter case, the 
small-r behaviour is reached in practice for r < 0.02 (Fig. 5); while for U/D ^> 1 Fig. 
7 shows that the critical A D r /U remains close to its Kondo asymptote of irr/8 over 
a wider r-range (up to r ~ 0.1), as one expects physically for U/D ^> 1 where charge 
fluctuations are relatively less significant. Note however that, even for U/D — > oo, the 
phase boundaries of the soft-gap Anderson and Kondo models coincide strictly only 
as r — > 0: the condition for mapping the former onto the latter via a Schrieffer- Wolff 
transformation — and hence the suppression of charge fluctuations in the Anderson 
model — is A X r /U < 1 (with A = min[D, U/2]), and from Eq. (2.10) is satisfied 
asymptotically on the phase boundary as r — > 0. Charge fluctuations in the Anderson 
model cannot therefore be neglected entirely even as U/D — > oo, save for r — > 0; as 
manifest in the fact (Fig. 7) that (AqD 7 " /U) c remains in general finite at the SC/LM 
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Figure 7. NRG phase boundary (AoD r /U) c versus r for U/D = 10 3 (solid circles), 
compared to that obtained via the LMA for U/D 3> 1 (solid line). Lines connecting 
NRG points are a guide to the eye only. Inset: results on an expanded scale, including 
the Kondo asymptote of wr/8 (dashed line). 



More broadly, Figs. 5 and 7 show the predicted LMA phase boundaries to be in very 
good quantitative accord with the NRG results for r < 0.3 or so; although with further 
increasing r the NRG boundaries increase more rapidly than their LMA counterparts. 
This brings us to the issue of how the phase boundaries behave as r — > |— (for r > |, 
both the NRG and LMA yield solely LM states for U > 0). Previous NRG studies for the 
soft-gap Anderson model [14,15] have suggested a divergent phase boundary as r — > \ — , 
although it has been argued in I that the data themselves do not reliably warrant such 
a conclusion. We have therefore re-investigated the matter more comprehensively, and 
conclude that the phase boundaries are indeed divergent as r — > |— . An incipient 
divergence is evident in the NRG results for U/D = 10~ 3 shown in Fig. 5, which has been 
extended up to r = 0.499; it is relatively weak, and fits well to logarithmic behaviour, 
(AqU t /U) c oc — m(| — r). Perhaps the most compelling evidence is however obtained 
by fixing r = \ precisely, and (for the chosen U/D and with D = 1 as the energy 
unit) progressively increasing the hybridization A : up to A values that are more than 
two orders of magnitude in excess of the critical A for r = 0.499, only LM states are 
observed and there is no hint of a SC state at r — \. Such behaviour is naturally not 
specific to U/D = 10~ 3 (although note that in Fig. 7 for U/D = 10 3 , the final r <\ 
NRG point is r = 0.49). 

The above behaviour as r — > |— is in contrast to that found within the LMA, where 
the phase boundary terminates at a finite value (although note ironically that a divergent 
phase boundary as r — > |— arises at the simple static mean- field level, see e.g. Fig. 9 of 
I). This discrepancy is plausibly explained by noting that while charge fluctuations are 
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partially included in the present (naturally approximate) LMA, it deliberately focuses 
on and successfully captures the strong coupling physics of the Kondo/spin-fluctuation 
regime that is asymptotically dominant for small-r. It does not therefore do justice to the 
charge fluctuations that, granted a divergent phase boundary, must become increasingly 
important as r — > \ — . Incorporation of such effects into the basic approach (beyond 
those already present) should however be possible, and is currently under investigation. 
We also add that a finite critical J c {j) as r — > |— has hitherto been reported in an 
NRG study of the symmetric soft-gap Kondo model [13]. This is now thought to be an 
artifact, however, and J c (r) is likewise believed to diverge as r — > \— [39]. 

Finally, we consider briefly the exponent for the vanishing of the Kondo scale as 
the SC/LM transition is approached from the SC phase, as reflected in the behaviour 
of the low-energy NRG scale u;*(r) (§3) whose critical behaviour may be expressed as 

w.(r) ex (l - P ^f) <r) ■ (4.1) 

Here, p = A X r /U (with A = min[D, U/2]), and the exponent x(r) is thus defined. 
The LMA result for x(r) is x(r) oc 1/r as r — > (Eq. (2.11)) and, while this can be 
established analytically only as r — ■> 0, numerical analysis shows the exponent 1/r to 
hold generally within the LMA (§6.1 of I). NRG results for x(r) are given in Table I; and 
as r tends to either or | we estimate a numerical uncertainty of 0(1) in the quoted 
values. (The data shown refer specifically to U/D = 10, although we do not believe 
them to be [7/D-dependent: the same values are obtained to within numerical accuracy 
for U/D = 0.1.) 

Table I. NRG results for the exponent x(r). 



x(r) 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.47 



19.6 10.8 7.6 6.2 5.43 5.14 5.2 5.85 8.55 12.2 



From Table I it is apparent that x(r) ~ 1/r as r ^ is indeed recovered. With 
increasing r, however, x(r) falls off less rapidly than 1/r, reaches a minimum at r ~ 0.3 
and increases again as r — > | where we suspect it diverges. While we cannot offer an 
explanation for this behaviour, save as r — > 0, we do not doubt at least its qualitative 
accuracy; and in view of the comments made above it is hardly surprising that the NRG 
and LMA exponents differ with increasing r. 



5. Single-particle dynamics 

We begin with a brief overview of spectral characteristics in both SC and LM phases 

for a representative r < |, and compare directly NRG and LMA determined spectra; 

r = 0.2 is chosen, together with a fixed ratio U/D = 0.1. Fig. 8a shows NRG results for 

i 

the single particle spectrum versus Co = u/A^ r , for four reduced interaction strengths 
U = U/Ap: U = 31.6 (LM phase), 19.0 (LM), 13.3 (SC phase) and 10.1 (SC); the 
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SC/LM transition occurs at U c ~ 13.7. Corresponding spectra obtained from the LMA 
for the same parameters are shown in Fig. 8b. 
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Figure 8. A£~ r D(u) versus u = w/A 1 " r . (a) NRG with fixed r = 0.2 and U/D = 0.1, 

and (from outside to inside) for: U = U/A^ = 31.6 (LM), 19.0 (LM), 13.3 (SC) and 
10.1 (SC); the critical U c = 13.7. (b) LMA spectra for the same parameters. 



With decreasing interaction strength in the LM phase the Hubbard satellites, 
centred on u = ±|£7, naturally move to lower energies; and in addition a narrow 
low-energy feature is seen to develop in the vicinity of the Fermi energy uj — 0. This 
becomes increasingly pronounced as U decreases toward the U c (r) at which the LM/SC 
transition occurs but, as seen from Fig. 8, D(u) vanishes as uj — > in the LM phase, 
with the asymptotic behaviour D(u) ~ \uo\ r (behaviour that is characteristic of the LM 
state for any r > 0). For U < U c (r) by contrast, the spectra exhibit the characteristic 
uj — > divergence D(uj) ~ M~ r that is symptomatic of the SC phase for any r < \. 
And with further decreasing interaction strength in the SC phase the Hubbard satellites 
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naturally become less distinct, and as U — > are subsumed into the central divergent 
resonance (see e.g. Fig. 11c of I). 

As is evident from Fig. 8, agreement between the NRG and LMA spectra is 
rather good (see also Fig. 15 of I); the principal difference being that the NRG 
Hubbard satellites in Fig. 8a are more diffuse than their LMA counterparts, and become 
increasingly so with increasing U. However this is likely to be an artifact of the 
NRG calculations. The logarithmic discretization employed in the NRG is designed 
to capture the low-energy physics; and on higher frequency scales appropriate to the 
Hubbard satellites, concomitant broadening of the discrete and relatively sparse NRG 
poles is well known to produce spectral overbroadening. We add moreover that the 
Hubbard satellites within the LMA already contain the effects of additional many body 
broadening beyond simple mean-field level, well known to exist for the normal r = 
Anderson model (see e.g. [20,32]), and shown in §7 of 1 to arise also for the soft-gap 
model; and on these grounds too the diffuse character of the NRG Hubbard satellites is 
unlikely to be physically correct. 

As discussed in §2 (see also I and [16]), the modified spectral function F(uj) = 
7tA [l +tan 2 (|r)] \u\ r D(u;), which removes the \oo\~ r divergence that is symptomatic 
of the SC state but unrenormalized by interaction effects, provides a much more acute 

probe of single-particle dynamics in the SC phase. This is illustrated in Fig. 9, again 

i 

for r = 0.2 and U/D = 0.1, where NRG results for J- (to) versus Co = cu/Aq - '' are shown 
for reduced interaction strengths: U = 8.0 (SC), 10.1 (SC), 13.3 (SC) and 14.2 (LM). 
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Figure 9. NRG spectra T{uj) = ttA [1 + tan 2 (§ r)\ \u\ r D{uo) versus Co = uo/A^ for 

fixed r = 0.2 and U/D = 0.1, and (from inside to outside) for: U = U / A^ = 8.0 
(SC), 10.1 (SC), 13.3 (SC) and 14.2 (LM). The SC/LM transition occurs at U c = 13.7. 

Note first from Fig. 9 that the spectral sum rule T{io = 0) = 1 (see §2), that 
generalizes to the SC state for arbitrary r < \ [16] the pinning condition familiar for 
the r = Anderson model, and encompasses the latter as a particular case, is well 



Anderson impurities in gapless hosts 



19 



satisfied in the NRG calculations; being obeyed typically to within an accuracy of 1-2% 
throughout the SC phase. 

As evident from Fig. 9, and in direct parallel to the normal Anderson model, 
J-(oo) in the SC phase contains both the expected Hubbard satellites and a generalized 
Kondo resonance that narrows progressively as U is increased towards the SC/LM phase 
boundary at U c (r); in fact the SC F(u) shown in Fig. 9 is qualitatively indistinguishable 
from that corresponding to the r = case on the 'low-resolution' scale shown (see e.g. 
Fig. 6 of [16]). For concreteness in what follows we define the Kondo scale ook = oo K {r) 
as the i-width at ^-height of ^(uo). The definition is of course somewhat arbitrary, the 
important point being that ook(j') should vanish as U — > U c (r)—, and within the NRG be 
proportional to the low-energy scale u;*(r) (with which it is hence in effect synonymous). 
Likewise, within the LMA, wk()") is proportional to the spin-flip scale u m (r) (§2), and 
thus correspondingly synonymous. 

In contrast to the normal Anderson model, the SC uj^(j-) vanishes at a finite U c (r), 
the effects of which are seen clearly in Fig. 9 where for U > C/ C (r)(~ 13.7) in the LM 
phase the Kondo resonance is absent, and as uo — > the LM !F(u) oc \uo\ r D(uo) vanishes 
with characteristic \uo\ 2r behaviour. The £7=13.3 (SC) and 14.2 (LM) examples included 
in Fig. 9, corresponding to 8 — \U /U c (r) — 1| ~ 0.03, also exemplify spectral evolution 
as the transition is approached from either side. They show, as argued within the LMA 
in I (see Fig. 12 therein), that the SC and LM spectra on either side of the SC/LM 
transition are near coincident on all energy scales save the lowest; and as 5 — > the 
LM/SC spectra coincide to arbitrarily low energies, in which sense the spectra evolve 
smoothly as the phase boundary is approached. 

6. Spectra: universal scaling 

Since the Kondo scale uj-^(r) - - and hence the width of the Kondo resonance in 
J- (to) - - vanishes as the SC/LM transition is approached from the SC side, one 
anticipates spectral scaling, with an r-dependent family of universal scaling spectra 
T = JF(cj/cjk(^)) obtained as U — > U c (r) — . This is indeed predicted by the LMA (§8 
of I), and is now investigated via the NRG. 

To this end we first consider the normal Anderson model, r = 0, where T{ui) = 
7iAoD(u). Here, spectral scaling has hitherto been examined via NRG for an asymmetric 
model [28], in the one-hole sector uo < 0; and for the symmetric model by a 
QMC/maximum-entropy study [30] at finite temperatures and for moderate interaction 
strengths U = U/A < 3n. NRG results for the T = symmetric Anderson model, 
J- (to) versus oo/ook, are shown in Figs. 10 and 11; they were obtained using A = 1.5 and 
Ag = 900, and universal scaling in terms of oo/ook is well established for U > 5ir. Fig. 
10 also compares the NRG scaling spectrum to that obtained previously [32] within the 
LMA. The inset shows the LMA results (Fig. 8 of [32]) arising when a particle-hole ladder 
sum is employed for the transverse spin polarization propagator that enters the LMA 
self-energies T, c (uj)\ while the main figure shows a simple modification detailed in [32] 
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Figure 10. For normal Anderson model (r = 0), NRG scaling spectrum J-{u>) — 
ttA D(uj) versus w/uk (dashed line), compared to that obtained via the LMA (solid 
line) as discussed in text. 

that merely eliminates the small spectral artifact apparent in the former at u/uk ~ 1-4, 
but otherwise produces no significant effect on the scaling spectrum. 

As seen from Fig. 10, the agreement between the LMA and NRG scaling spectra for 
r = is rather remarkable over the full frequency range. In addition to spectral pinning 
at the Fermi level (u — 0), the LMA successfully captures both the characteristic low- 
frequency Fermi liquid behaviour D{u) — D(0) ~ (uj/ujk) 2 and, at larger uj/uk, the 
Doniach-Sunjic (DS) law [27,28,30] indicative of the orthogonality catastrophe, whereby 
D(uj) ~ (\uo\/uo K )^ a with a — 1 — 2(5 /tc) 2 and 5 = ix/2 the Fermi level phase shift; ie 
D{yj) ~ (|cj|/c<jk)~2 . The existence and, indeed, dominance of the latter is apparent in 
Fig. 11, where a DS fit of form 

= a + b(\Lu\/u K )-^ (6.1) 

is made to the NRG data (a fit to the LMA is essentially equivalent). As seen from the 
inset to Fig. 11, DS behaviour persists out to uj/uj k = 30 (and naturally well beyond), no 
deviation from the form Eq. (6.1) being apparent; while the main part of Fig. 11 shows 
clear DS tail behaviour to extend down to u/uk — 0.5. By contrast, the low-frequency 
Fermi liquid behaviour persists only up to uj/uj k ~ 0.2, where a rapid crossover to DS 
behaviour begins; this too is seen in Fig. 11 where extrapolation of the Fermi liquid 
portion of the spectrum (D(u) — D(0) oc [c<j/c<jk] 2 ) is made, producing a Lorentzian 
J~{<jS) that, save for uj/uj-k < 0.2, fails entirely to capture the scaling spectrum. 

We turn now to r > 0, to consider scaling of the SC as the SC/LM transition 

is approached. Specific NRG results shown below are obtained for fixed U/D = 0.1 

i 

(as in §5); but we note that while the particular value of U c (r)(— (U/ A^ r ) c ) at which 
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Figure 11. For r = Anderson model, NRG T(oj) = ttA D(uj) versus w/wk including 
fits to the DS tails (Eq. (6.1), dashed line) and to the Fermi liquid portion of the 
spectrum (D(uj) — D(0) oc (w/wk) 2 , dotted line). The inset shows the continued 
persistence of the DS fit. 

the transition occurs naturally depends upon U/D (see e.g. Fig. 6), variation of this 
parameter produces no detectable influence on the scaling spectra themselves. For 
r = 0.2, Fig. 12 shows the resultant NRG F(w) for U = 8.0, 10.1 and 13.3 (as in 

Fig. 9); the critical U c (r) ~ 13.7. The left-hand inset to the Fig. shows the central 

i 

portion of the Kondo resonance, versus Co = uj/A^ r on an 'absolute' scale, to illustrate 
its rapid narrowing with increasing U, and ultimate vanishing as U — > U c — ; the same 
behaviour is clearly evident within the LMA in Fig. 18 (inset) of I. By contrast, the main 
part of Fig. 12 shows the scaled F{oo) versus w/wk, from which universality is indeed 
apparent; the scaled spectra for U =10.1 and 13.3 in particular are to all intents and 
purposes indistinguishable. Note also that while U c (r) is finite, the Hubbard satellites 
are naturally projected out of the scaling spectrum since oj-^ij-) vanishes at the phase 
boundary. 

The LMA predicts two essential characteristics of the SC spectra (§8 of I). First, 
that for r > the leading u/u^r) ~~ * behaviour of — 1(= F(uo) — ^"(0)) is 

T{uj) — 1 oc [IcjI/cjk] 7 ^ with 7(r) = 1 — r (and a proportionality constant that naturally 
vanishes as r — > 0, see §8 of I). Cusp behaviour is indeed seen in the r = 0.2 NRG 
results of Fig. 12, and numerical analysis confirms an exponent of j(r) = 0.8 (asymptotic 
behaviour that in practice sets in for |u;|/c<jk < 0.01). Second, that the wings of the 
scaling spectrum should contain generalized DS tails, of form F(uo) ~ [IcjI/cJxfV)] ^ 
with an exponent u(r) — \ — r. This is likewise found in the NRG calculations, as seen 
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Figure 12. For r = 0.2, NRG scaling spectrum = 7rA [l + tan 2 (f r)] \u\ r D(uj) 

versus uj/uj k ; with fixed U/D = 0.1 and for U = U/A^ 7 = 8.0, 10.1 and 13.3 (the 

i 

critical U c (r) = 13.7). Left inset: corresponding data versus u — uj/A^ t , to illustrate 
vanishing of Kondo scale as U — > U c (r)—. Right inset: DS fit to the scaling spectrum 
(dashed line), Eq. (6.2) with tail exponent v(r) = \—r. 



from the right-hand inset to Fig. 12 for r = 0.2, where a DS fit of form 

= a + b[\u\/u K (r)]-^ (6.2) 

is made. The fit is excellent (extending down to |w|/c<jk ~ 0.5) and, as for the r = 
model, the DS tails are again seen to dominate the scaling spectrum. The above 
behaviour is not of course confined to r = 0.2. Fig. 13a shows NRG results for T{uS) 
versus uo/oo^r), for r = 0.1, 0.2 and 0.3 (from bottom to top); the inset shows the 
same data on a reduced scale \u\/u-^{r) < 1. Corresponding LMA results are shown in 
Fig. 13b (obtained using the particle hole ladder sum in and with no attempt 

to avoid the spectral anomaly apparent around |cj|/wk ~ 1~2). In all these cases, 
for both NRG and LMA spectra, the cusp exponent is found to be 7(7") = 1 — r; 
and DS tails of form Eq. (6.2) again arise. We also add that the DS exponent of 
v{r) — \ — r reflects the fact that we consider T{uS) oc \uj\ r D{uj) which, unlike D{uj) 
itself, is a universal function of u/uo^r) as the transition is approached from the SC 
phase. However [uo^{r)] r D{uo) oc {\uo\ / uo-^{r))^ r J 1 {uo) does exhibit scaling, and we note 
that its DS exponent is thus u(r) + r — \ — precisely as for the normal Anderson model, 
and again symptomatic of the fact that the latter consitutes a particular case of more 
generic SC behaviour. 

In qualitative terms, both NRG and LMA results for the universal T{w) in Fig. 
13 show narrowing of the central portion of the Kondo resonance with increasing r 
(insets to Fig. 13); together with a relative 'flattening' {y{r) — \ — r) and increase in 
intensity of the DS tails. Quantitatively, the agreement between NRG and LMA for 




Figure 13. Scaling curves = ttA [1 + tan 2 (|r)] \uj\ r D{lj) versus w/w K (r) for 

(from bottom to top) r = 0.1, 0.2 and 0.3; as obtained via (a) NRG, and (b) LMA. 
Insets: same data on a reduced scale; from outside to inside the curves correspond to 
r = 0.1, 0.2 and 0.3. 



J-{oj) becomes poorer with increasing r — not surprisingly in view of the discussion of 
the phase boundaries given in §4, and reflecting the fact that the specific LMA we have 
considered focuses on the spin-fluctuation physics which is asymptotically dominant as 
r — > 0. Explicit comparison between NRG and LMA results for the universal is 
made in Fig. 14 for r = 0.1; the most obvious difference being the relative offset in the 
DS tails (each of which fits very well to the form Eq. (6.2)), which increases with r (Fig. 
13). This shows up in significant differences between the NRG and LMA a-coefficients 
in Eq. (6.2); although the corresponding ^-coefficients are much closer, reflected in the 
fact that the DS tails in Fig. 14 are near parallel. 
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Figure 14. Direct comparison of r = 0.1 scaling curves ^F(u>) = 
ttA [l + tan 2 (fr)] \uj\ r D(lu) versus uj/lj k : NRG (upper) and LMA (lower). 

7. Summary 

In this paper we have studied the symmetric soft-gap Anderson impurity model [9,14- 
16,31,32]. The model exhibits a rich spectrum of physical behaviour on account of the 
underlying quantum phase transition between strong coupling and local moment states, 
and contains the 'normal' Anderson model as a particular limiting case. Our primary 
aim has been to make detailed comparison between two non-perturbative approaches: 
the NRG method, which provides essentially exact numerical results for the problem 
[14,15], and a recently developed many-body theory, the LMA [31,32]. 

We find excellent agreement between the two approaches, both for static properties 
such as the ground state phase diagram, as well as for dynamical properties embodied 
in zero temperature single-particle excitation spectra. In particular, specific predictions 
arising from the LMA [31] are found to be well borne out by NRG results. These 
include scaling characteristics and asymptotic behaviour of the phase boundaries, Fermi 
level pinning of the modified spectral functions throughout the SC phase, and universal 
scaling of single-particle spectra upon approach to the SC— >LM transition where the 
Kondo scale characteristic of the SC phase vanishes. Differences between the two 
approaches arise only as r — > |; they may be ascribed to an incomplete treatment 
of charge fluctuations in the specific LMA we have considered in practice, which seeks 
primarily to capture the spin-flip physics that dominates the regime of well-developed 
local moments (whether in the SC or LM phases). Quantitative agreement between 
NRG and LMA results is however remarkably good in general, the universal scaling 
spectrum of the normal Anderson model (Fig. 10) providing a particular example. The 
results obtained here provide in our view strong encouragement to further study of the 
present — and related — quantum impurity models, via both the NRG and LMA, and 
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including consideration of finite-temperature spectral and hence transport properties. 
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